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Abstract. Correlations in photodetection signals from quantum light sources are 
conventionally calculated by application of the source master equation and the 
quantum regression theorem. In this article we show how the conditioned dynamics, 
associated with the quantum theory of measurements, allows calculations and offers 
interpretations of the behaviour of the same quantities. Our theory is illustrated for 
photon counting and field-amplitude measurements, and we show, in particular, how 
transient correlations between field-amplitude measurements and later photon counting 
events can be accounted for by a recently developed theory of past quantum states of 
a monitored quantum system. 


1. Introduction 

The experimental observation in 1956 [I] by Hanbury-Brown and Twiss of temporal 
correlations in photo detection signals constitutes a defining moment for the field of 
quantum optics. While the signals measured were compatible with classical field 
amplitude fluctuations, in a matter of few years it became evident that temporal 
correlations in signals from quantum light sources, such as a single atom, would 
defy classical interpretation. Many scientists have since then contributed to these 
developments and studies of how correlations in photodetection records reveal the 
quantum properties of the emitter, of the light field, and even of the measurement 
process itself. Today, optical detection still holds a prominent role in both fundamental 
tests of quantum theory and in efforts to explore quantum phenomena, e.g., in precision 
sensing and secure communication protocols. 

Correlation functions play an important role in the analysis of time dependent 
processes with applications ranging from demographics and finance to engineering and 
physics. Stochastic realizations of different signals can be characterized by, e.g., their 
mean values and variances while their temporal dynamics is captured by correlations 
between their values at different times. One may thus evaluate products of two signals 
at different times, nii(ti)m 2 (t 2 ), and define the temporal (two-time) correlation function 
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of the signals as the average of this quantity, 777.1(^)7712(^2)) over different realizations 
of the process. While the expression deals with classical signals mi(t) and m 2 (t), they 
may have their origin in measurements on a quantum system, and hence we must apply 
quantum theory to analyze and make predictions for the signal correlations. 

A central position in this research is held by Glauber’s photodetection theory M, 
which explicitly recognized that in the process of light detection, the measurement signal 
is a fluctuating classical current, i(t), formed by electrons that have been excited by the 
absorption of photons. The mean value i{t) is proportional to the quantum expectation 
value (a){t)a{t)), but to determine the correlation function for t\ < t 2 , we 

must incorporate the effect of subsequent individual photon-absorption events, leading 
to the normal- and time-ordered quantum expectation value (dt (t^a} (t 2 )d(t 2 )d(ti)) of 
the field creation and annihilation operators. 

Glauber’s photodetection theory is formulated in the Heisenberg “operator picture” 
of quantum mechanics, and the evaluation of correlation functions traditionally applies 
the master equation and the quantum regression theorem for the evaluation of the 
observables of the light emitting system [ 5115 ] . In this article we show how photodetection 
theory can also be naturally formulated within the general quantum theory of 
measurements, which evaluates the the joint probability distribution P(? 77 i, m2, £2), 

that photodetection signals acquire the values 77 i 1 (t 1 ) = mi and 7772(^2) = tt 7 2 , by careful 
evaluation of the back action on the emitter state due to the first measurement and 
its subsequent evolution until the second one. When one knows the joint probability 
distribution, the correlation function follows by 

7771(^1)7772(^2) = {m 1 m 2 )P(mx,t 1 ; m 2 ,t 2 ), ( 1 ) 

mi,m2 

where sums can be interchanged by integrals in the case of continuous spectra of 
measurement outcomes. 

The quantum measurement theory approach accounts explicitly for the back action 
on the state due to measurements on the system in a manner similar to the stochastic 
quantum trajectory description [ 7 ], Thus, it both reproduces the results of the quantum 
regression theorem and sheds light on the origin of the temporal correlations in optical 
detection records. We demonstrate here that the recent past quantum state formalism [ 8 ] 
provides a similar back action, but on the state prior to a measurement event and thus 
offers a “symmetric intuition” to the inherently asymmetric case of cross correlations 
between intensity and amplitude measurements. 

The article is organized as follows: I11 section [ 2 ] we introduce the master equation 
and recall how the quantum regression theorem enables the calculation of photodetection 
correlation functions. In section [ 3 ] we review some aspects of the general quantum 
theory of measurements that allows an alternative formulation of correlation function 
calculations, which we extent to the past quantum state formalism in section [fl In 
section [ 5 ] we use this alternative formulation to determine and explain interesting 
correlation functions between intensity and amplitude measurements of light emitted 
by two- and three-level atoms, and section [ 6 ] concludes the article. 
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2. Correlation functions in photodetection theory 

2.1. The master equation and the quantum regression theorem 

A light-emitting quantum system can be described by a density matrix p which, under 
the assumption of the Born-Markov approximation, i.e., in the case of weak coupling 
to memoryless reservoir modes, obeys a linear master equation, 

P = £(h (2) 

where 

Cp =2-[if,p] + Z Cr.pCi - i(C*C„P + pC' n C„). (3) 

H is the Hamiltonian of the system (possibly driven by time-dependent external classical 
fields) while the operators C n account for dissipative couplings to the environment of 
the system. 

If the density matrix elements pij are arranged in a vector, £ can be represented 
as a matrix, and if this matrix has constant coefficients, equation ([ 2 ]) is formally solved 
by matrix exponentiation. The propagator of the master equation is a little more 
complicated than a simple matrix exponential if the equation is time dependent, but we 
shall nonetheless denote it by the symbol e Cr and use the expression 

p(t + r) = e Cr [p(t)] (4) 

for the solution to equation (J2]) . 

As stated in the Introduction, Glauber’s photodetection theory reflects how a 
detection signal is obtained by absorption processes in the detector. The average product 
of intensities and amplitudes measured at different times is thus formally related to 
the expectation value of a normal ordered product of held creation and annihilation 
operators, a^(t) and a(t'). The mean intensity at time t is given by (cd(t)d(f)), while 
the intensity-intensity correlation function, i.e., the expectation value of the product 
of intensity measurements at two times £ < t 2 is G^(ti : t 2 ) = (ti)a) (t 2 )a(t 2 )a(ti)), 
where the superscript ( 2 ) refers to the quadratic dependence on intensity. 

Like in classical electrodynamics, the electric and magnetic fields can be expressed 
as functions of the charge and current distributions, and in [9] it is described in detail 
how the held creation and annihilation operators develop components proportional to 
the raising and lowering dipole operators of the light emitting system. For a transition 
in an atom between a definite pair of excited and ground states \e) and | g), the atomic 
emission of photons with rate 7 is represented by the operator C = ^\g)(e\ in (J3]) , 
and hence the relevant held annihilation and creation operators, a and a\ are replaced 
by the atomic transition operators a = \g)(e\ and cfl = \e)(g\. The two-time intensity- 
intensity correlation function is thus proportional to the correlation function of atomic 
operators, 


G (2 \t 1 ,t 2 ) oc (d t (ti)a\t 2 )a(t 2 )a(ti)). 


(5) 
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The two-time expectation values in fl5J) is on the form (0i(t)0 2 (t + r)0 3 (t)) = 
Tr(0i(f)0 2 (£ + r)0 3 (f)p(t)) = Tr(0 2 (t + r)0 3 (t) p(t)Oi(t)), where the central operator 
0 2 (t + r) (= <xl(t 2 )cx(t 2 )) can be expanded on the complete set of dyadic operators 
(|j)(i|)(t+r). The expectation values of these operators are the density matrix elements, 
Pij = ((|j)(i|)), which obey the linear set of coupled equations defined by the master 
equation (J2]). The same is, however, not true for the operators (|j)(i|) themselves, 
and noise operator terms must be added to the equations to ensure the preservation of 
commutation and uncertainty relations |5],[6j. In the Markov approximation, however, 
these noise operators are uncorrelated with all system observables at earlier times, 
and this implies that two- and multi-time operator correlation functions of the form, 
(Oi(t)(\j)(i\)(t + r)0 3 (t)), indeed, evolve with the time argument r according to the 
same linear set of coupled equations as ((|j)(i|)(t + r)), i.e., as pij(t + r), as long as 
r > 0 and Oi(t ) and 0 3 (t ) represent operators, evaluated at earlier times. This is 
the Quantum Regression Theorem png, and it offers an effective means to determine 
correlation functions of operators for Markovian open quantum system. 

For comparison with experiments, we are interested in steady state dynamics, and 
hence only in the dependence of the correlation functions on the time difference r. We 
thus assume that p[t ) = p st , the steady state solution of the master equation, and, for 
r > 0, the quantum regression theorem yields the formal solution mm 

<0i(t)0 2 (* + r)0 3 (t)) = Tr (d 2 e CT [0 3 p st d 1 ]) . (6) 

The master equation is linear and preserves the trace and hence, for any matrix /x, 
e £r [/i] —> Tr (p)p st for r —> oo. This implies that the correlation function (J6]) regresses 
from the steady state mean value of the operator product (0i0 2 0 3 ) st at r = 0 to the 
product of steady state mean values (0 2 ) st (0i0 3 ) st . In quantum optics it is convenient 
to normalize the correlation functions by this product of the steady state mean values. 
The resulting reduced correlation functions, denoted by a lower case letter g, then 
approach unity for large r. 

2.2. Intensity-intensity and intensity-amplitude correlations 

With the choice of operators 0\ = cr\0 3 = a, and 0 2 = (fid, the intensity-intensity 
correlation function G^(r) (J5]) is given by ([6]). With these operators, 0 3 p st 0\ = 
pf e \g){g\, i.e., the atomic ground state density operator, multiplied by the steady state 
excitation probability pf e . In this case, the evolution by exp(£r) merely solves the 
transient evolution of the atomic density matrix, starting from the ground state, and 
since 0 2 = |e)(e| is the projection operator on the excited state, the reduced correlation 
function evaluates to g (2) (r) = pf e p e e(r)\ g /(pf e ) 2 = Pee{r)\ g /pf e , where p ee (r )\ g denotes 
the excited state population of a system, evolved from the ground state at r = 0. This 
result is readily understood from the quantum trajectory picture [711111112] : The photon 
detection at time t is accompanied by an atomic quantum jump into the ground state, 
and the probability to register another photon from the same atom at the later time 
t + r is proportional to the excited state population conditioned on the first detection. 
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Figure 1 . A setup for joint photon counting and homodyne detection. Light 
transmitted by the beam splitter along the lower path is detected as discrete photon 
counting events, while the light reflected along the upper path is mixed with a strong 
local oscillator and the difference between the two counting signals represents the 
continuous homodyne monitoring of the emitted field amplitude. Averaging the noisy, 
time dependent homodyne signal in time windows around all counting events in the 
detector Di, we obtain the amplitude-intensity correlation function. 


The particle and the wave description of optical phenomena competed as the 
dominant descriptions of light until Planck and Einstein introduced the dual quantum 
nature of light more than a century ago. In beautiful, more recent demonstrations of 
the joint wave and particle properties of light [TTHH4] . the radiation emitted by an atom 
was split in two components, see figure [D for which the photon number and the field 
amplitude were detected, respectively. These experiments revealed strong temporal 
correlations between the measurement outcomes. For a positive time delay, r > 0, 
between the counting events and the amplitude measurement, such amplitude-intensity 
correlations are captured by the time- and normal-ordered coherence function m 


<7 {1 ' 5) (t>0) 


+ r) + e lip a{t + r))cr(t)) 
Pte(^P% + z-' lip P%) 


(7) 


where (p is the phase component measured by the homodyne detector, and where the 
superscript (1.5) refers to the “intensity 1 " 5 ”-dependence of the function. The numerator 
in this expression is represented by (|6j), with 0\ = di, 0 3 = a, and 0 2 = + e _1</, d, 

and its behaviour can be understood from the quantum jump of the steady state density 
matrix into the ground state at time t, and the evaluation of the expectation value of 
the amplitude operator 0 2 during the subsequent transient evolution of the conditional 
state. 

In the experiments [TTIlfld] . the homodyne amplitude signal shows correlations also 
with later photon counting events. The measured correlations involve the same operators 
Oi as before, but for r < 0, the time- and normal-ordering procedure leads instead to 
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Figure 2. Amplitude-intensity correlation function 0 and © for a resonantly driven 
two-level atom determined by the quantum regression theorem (blue solid line) and 
by the past quantum state (red dashed line) (THU) . The parameters are: = 5q and 
ip = 7 t / 2 . 


the expression [15 ] 


g (L5) (r < 0) 


- r)a(t - r)) + e - r)a(t - r)a(t)) 

pte(^pfe + ^P%) 


In the evaluation of (JHJ), the time evolution operator in fl6j) is applied to products of 
matrices and ap st that do not permit separate interpretation as quantum states, 
and unlike the case of r > 0, we do not have a simple interpretation of the expression 
in dHJ) in terms of the transient conditioned dynamics of the system. Wiseman has 
shown |l6j, however, that the weak value formalism DU for the post selected average 
of a weakly perturbing measurement, indeed, accounts for field-amplitude correlations 
with later count events. That formalism is related to the recently developed theory 
of past quantum states [8], and as described below, the general quantum theory of 
measurements and the past quantum state, allow both calculation (the red dashed curve 
in figure [2]) and interpretation of how field-amplitude measurements are correlated with 
later detector click event. 

In figure 0 we show the value of the normalized correlation function g^' 5 \r) for a 
two-level atom for both positive and negative r, and it is interesting that despite the very 
different formal expressions, the evaluation of (JTJ) and (jHJ) yields a symmetric correlation 
function around r = 0. Since intensity-intensity correlation functions concern fully 
equivalent detection events, they have to be symmetric, but as we shall see below, 
intensity-amplitude correlation functions may, indeed, be very asymmetric. 
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3. Quantum measurement theory and field correlations 

Equation (J2]) describes the average behavior of an unobserved light emitting system, 
while observation of the emitted radiation will yield a conditioned dynamics, described 
by the general theory of measurements. Formally, any observation of a quantum system 
is described by a positive operator valued measure (POVM), i.e., a set of operators fl m , 
with Y.m&lrP'rn = 7, where I is the identity matrix pJlIIHJ- The operators Q m account 
for the joint evolution of the system and an appropriate meter system and read-out of 
the meter in the states that we associate with the outcome results m. The probability 
to get a definite outcome m is given by the expression, 

P(m) = Tr(pO^Q m ), (9) 

and the (normalized) state of the system, conditioned on this measurement outcome, is 

Pm = n m pCl\n/P{m). (10) 

If {fl m } constitutes a complete set of orthogonal projection operators, the POVM 
formalism yields the conventional Born rule and von Neumann projection postulate. 

The POVM formalism also gives access to the probability that a measurement, 
described by POVM operators Ct m yields mi at time p and a subsequent measurement 
at t 2 on the same system, described by the same or a different set of POVM operators, 
yields the outcome m 2 . Such joint probabilities follow Bayes’ rule, and can be written 
as 

P(mi,ti,m 2 ,t 2 ) = P(m u ti) ■ P(m 2 ,t 2 \7n 1 ,t 1 ). (11) 

Assuming steady state conditions, the unconditional probability P(mi) =Tr(O mi p st Oj ni ) 
for the first event to occur does not depend on the time argument t\, and the state of 
the system, conditioned on this outcome is p mi (ti) = O mi p st O) ni /P(mi). Propagating 
the density matrix further, until t 2 , by the usual master equation with p mi (t\ ) as the 
initial condition, we find 

PmAh) = e £(t2_tl Vmi(P)- (12) 

For the joint outcome probability distribution, we thus obtain 

P(m 1 ,t 1 -m 2 ,t 2 ) = r &(H m3 e /:(ta-tl) [n rni p pt n| ni ]nJ n2 ), (13) 

where acts on the matrix object contained within the square brackets. 

Equation (U3l) provides the outcome joint probabilities and offers a microscopic, 
dynamic view of a conditionally evolved quantum state. This differs from ((6j) which 
directly provides the average correlation function of the physical observables, e.g., of 
the field and intensity operators. They are connected, however, by using equation (fT3l) 
for the calculation of weighted mean values of the outcomes: If the labels mi and jn 2 
denote the actual values of the physical observables measured at times t\ and t 2 , the 
operator correlation function of interest acquires the form (1), 

m 1 m 2 = 5Z (m 1 m 2 )P(m 1 ,ti,m 2 ,t 2 ). 

mi,m 2 


( 14 ) 



3.1. Photon counting and homodyne detection 


Restricting ourselves to the case of detection of light emitted from a quantum system, 
photon counting associates two POVM operators to the click and no-click detector events 
during an infinitesimal time interval di. The observation of a photon from a two-level 
atom is described by IR = y/'ydta, while the absence of a detector click is described, to 
first order in d£, by Cl 0 — (I — ^cFcr), reproducing the respective probabilities Pee'ydt 
and 1 — Pee'ydt as well as the final atomic states after the detection events. The atomic 
states are notably not orthogonal for the two held detection outcomes. The average 
behavior of the system, p —> subject to the two outcome possibilities is 

exactly the one described by the damping terms in the master equation (J2D with a 
single C = y/^cr. 

In homodyne detection, the light signal is mixed with a strong classical held of phase 
</?, and the intensity of the interfering helds is recorded by two photon counters, as shown 
in figure [D The difference between their counting signals is represented by a continuous 
quantity x, which in the absence of the emitter is represented by a Gaussian distribution, 


Po(x) = 


yj 2ir/dt 


exp(— x 2 dt/2), with a suitable normalization of the outcome argument 


x p2, 16j. In the presence of the two-state emitter, the outcome is governed by the 
family of POVM operators, 





exp 



(1 + axl ty/^e — 


(15) 


To hrst order in d£, one readily verihes the POVM property / Of O^da; = / and that the 
average change of the state, p —> f Vt x pVt\dx, reproduces the effect of the damping term 
C = yfyd in (|3J). 

The mean intensity measured as the average number of detection events per time, 
n/dt , has an expectation value given by the weighted mean of the possible outcomes of 
the photon counting measurement, ii = X) n =o l n-P(n) = X) n =o l w-Tr (ClnpClP) = pf^dt, 
where the last expression applies for the two-level atom in steady state prior to the 
measurement. Similarly, in homodyne detection, the one-time mean amplitude of the 
held is proportional to (x) — f x ■ P(x)dx = f x ■ Tr(Q x pfil)dx = y/j Tr(e~ lv, ap + 
e*pfi f ) = y/7(e- i *P?g + e* > p$ l ). 

The two-time intensity-intensity correlation function is proportional to 

712 = 0,1 nin 2 P(ni, tp n 2 , t 2 ) = P(n i = l,ti;n 2 = 1,^) with both of the POVM oper¬ 
ators in equation (TT3l) pertaining to photon counting. Applying the expressions for a 
two-level emitter, we thus obtain P(n\ = l,ti;n 2 = l,t 2 = ti + r) = ('fdt) 2 pf e p ee (r)\ g , 
which is proportional to G^(r > 0), identified in the previous section. In exactly the 
same way, we can determine the probability for the joint detection of a photon at time 
ti and the subsequent amplitude outcome x at time f 2 , P(n, t\, x, t 2 ), and we can deter¬ 
mine the average J2 n f nxP(n, tp, x,t 2 )dx. For positive time lag r — t 2 — ti this yields 
Pf e ( e ~ lip Peg{x)\ g + e lip p ge (r) \ g )^y/ldt which is indeed proportional to G^- 5) (r) Q. 
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4. Amplitude—intensity correlations and past quantum state 

We shall now address the joint probability of a homodyne detector outcome x at the 
earlier time t\ followed by the counting of a photon at the later time t 2 , i.e., we take 
hlmi and Cl m 2 kin in (1131) , 

P(x,ti,n, t 2 ) = P(x,t i) • P(n,t 2 \x,t 1 ) 

= Tv(n n e c ^\Q xP st Qlp n ). (16) 

To evaluate the amplitude-intensity correlation function, we multiply the probability 
by n and x, and we sum/integrate over their possible values. This will retain only the 
n = 1 component and the contributions that are first order in x in the expressions for 
QW- Thus, even though the Q x operators differ only infinitesimally from the identity 
and the back action of the x measurement is weak, when we calculate the integral over 
x, only the terms oc dt^ye~ lip a and dt^fe^a^, cf. (15), contribute to the correlation 
function average over the very nosy homodyne signal. 

Rather than merely reproducing existing results with a different method, the 
purpose of this work is to analyze the problem and gain new insight into the results. 
We thus rewrite equation (TT6j) for the joint probability of the outcomes, as P(n,t 2 ) • 
P(x, ti\n, t 2 ) to emphasize the conditioning on the later photon counting rather than 
the earlier homodyne detection event. Indeed, such retrodictive conditioning can be 
applied to any type of measurements, cf. the past quantum state formalism [8j , applied 
in [T91I20]. We shall briefly recall this formalism before specializing to the calculation 
and interpretation of amplitude-intensity correlations. 

The density matrix evolves according to a completely positive map, and hence, 
for any operator /i, we can formally write the action of the propagator of the master 
equation as a Kraus map PHttEJ, e £( * 2 tl} [/i] = Kq/iKI, where E Q K} y K a = /. 
The operators K a represent the unobserved time evolution between t\ and t 2 , and they 
depend of course on the explicit form of the master equation and the propagation time, 
but for our purpose we shall not need their explicit form. The cyclic properties of the 
trace allows us to transform equation (fT3l) into 

P(mi, ti; m 2 , t 2 ) = ^ K l^ m prm K a )- (17) 

a 

It is now natural to define the operator product E{t\) = E a K a . The order 

of the operators K' a and K a has been swapped compared to their action on the time 
dependent density matrix, and the continuous equation of evolution for E(t) is similarly 
obtained by swapping the Lindblad operator terms (J3J) in the master equation (J2J), 

(ipi = '-pH, S] + £ ClEC n - \(ClC n E + EClC n ). (18) 

We solve the equation for E(t) backwards in time from t 2 to t\ with the boundary 
condition, E(t 2 ) = 

The correlations in the joint probability distribution imply that knowledge of the 
later measurement outcomes alters the conditional probability for the earlier ones. We 
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can formally express this by P(mi, tp, m 2 , t 2 ) = P{jn 2 ,f 2 ) • P{jn\, £i|m 2 , i 2 ). We now 
determine the first factor -P(m 2 ,t 2 ) by summation over the index mi in (ITT)) and the 
conditional probability then acquires the form 


P(m 1 ,ti\m 2 ,t 2 ) 


Tr(a» 1 p(f 1 )nj„ 1 g(ti)) 
Em Tr(Q m p(t 1 )Q! n E(t 1 ))’ 


(19) 


where the matrix E(t\) and the density matrix p{t\) play symmetric roles, representing 
how the outcome probabilities are correlated with the prior and the posterior dynamics 
and measurements of the system, respectively [ 8 ]. We are interested in steady state 
correlation functions, and since we eventually want to account for only the two-time 
correlation with a single later photon detection event, p{t\) = p st . When we solve 
equation (fT 8 |) . we obtain predictions in full agreement with the quantum regression 
theorem results. 

As E{t) solves an equation which is similar to the master equation for p(t), we 
can exploit our intuition for density matrix evolution and infer that for the two-level 
atom, E(t) performs an evolution backward in time, starting from the excited state, 
which is similar to damped Rabi oscillations of the conventional density matrix. This, 
indeed, gives a qualitative explanation of the red dashed curve in figure El The evolution 
of E(t ) involves population loss out of the excited state by the Lindblad damping 
term — ^(C^CE + EC^C) with C = but also population feeding into the excited 

state by C^EC. Unlike the master equation (J2]) , equation (TT 8 |) does not preserve the 
trace, and unlike the density matrix converging towards its steady state value, p st , E(t) 
converges towards the identity matrix for long (negative) propagation times. It is thus 
not evident that the correlation function in figure El should be completely symmetric in 
the time argument, but it follows by a closer inspection (and solution) of the equations 
of evolution. 


5. Amplitude—intensity correlations from three-level atoms 

We shall in this section study the fluorescence from optical transitions in three-level 
atoms. These systems have alrady been analyzed PH331EDI22], by the quantum 
regression theorem, and they show a number of features, that we can explain with 
our theory. 

5.1. Ladder system 

Let us consider a three-level ladder atom with energy levels as drawn in figure [3](a). Its 
Hamiltonian can be written as 

H = — — (flidoi + LI 2 CT 32 + H.c.), (20) 

where a t] = \i){j\, ( i,j = 1 — 3), fli )2 are the (real) laser Rabi frequencies. Quantum 
jump operators of Q are CT = y / qTdi 2 and C 2 = with decay rates yi )2 , and 
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Figure 3. A three-level atom driven by two resonant laser fields. In the ladder 
configuration (a), the fields are of comparable strength while the upper state |3) has 
longer lifetime than the intermediate state |2). In the Id-configuration (b), the driving 
field and the decay rate on the |2) o |3) transition are much stronger than those on 
the |1) O |3) transition. 



Figure 4. (Color online) Amplitude-intensity correlation on the lower transition of 
a single three-level ladder atom by invoking the quantum regression theorem (blue 
solid line) and past quantum state (red dashed line). The parameters are: 121 = 371, 
1^2 = 571, 72 = O.271, 7 P h = 5 x 10 _5 7 i , and ip = 7r/2. 


C 3 = 1 /7ph(<3'33 — 022 — 071 ) with a dephasing rate 7 P h of the upper state with respect 
to the two lower states. 

We have calculated the correlation function go 1 ’ 5 ' between the intensity and 
amplitude signals in the frequency range of both the lower and the upper transition 
in the system. The quantum regression theorem, the POVM formulation, and the past 
quantum state formalism naturally all give the same quantitative predictions for this 
correlation function, shown in figure |4j 
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Unlike the two-level atom, the three-level atomic system shows a striking temporal 
asymmetry of the correlation function [15] . In particular, for the lower transition the 
frequency of the correlation function oscillations for the negative time delay is twice the 
one for positive time delay (see figure EJ) . The past quantum state reasoning explains 
why this happens. For r > 0, the photon counting event makes the atom density 
matrix evolve forwards in time from the ground state |1) at time t. For r < 0, the held 
amplitude can be calculated from the matrix E, evolving backwards in time from the 
excited state of the monitored transition, i.e., the intermediate state |2) in the ladder 
configuration. In equations ([3]) and (fT8|i the evolution due to the Hamiltonian is the 
same (modulo a sign), and if we disregard the dissipative terms, the evolution is readily 
analyzed in terms of the Hamiltonian eigenstates, 

1+) = ^(cos0|l) - |2) + sin6»|3)), (21) 

|0) = — sin 6*| 1) + 008 6*13), (22) 

h) = ^=(cos0|l) + |2) + sin0|3)), (23) 

with corresponding eigenenergies 

£+ = t^Hr, £o = 0, S- — ——Hr, (24) 

where Hr = \JVt\ + fi| and tan0 = 

After a photon counting event, the conditioned evolution starts with the atom in 
state 11), which is a superposition of all three dressed eigenstates. It subsequently evolves 
with relative phase factors that oscillate with the spectral separation Hr/ 2, which then 
yields the oscillation frequency of the density matrix elements. The backward evolution 
of E(t), on the other hand, starts from the state |2), which can be expanded as a 
superposition of only the |+) and |—) eigenstates which are separated by Hr. The 
matrix elements of E(t) thus evolve with a relative frequency of Hr, which therefore 
governs the twice faster oscillations of the two-time correlation function for negative r. 

5.2. V-system 

The Hamiltonian of the V -configuration depicted in figure [3](b) reads 

H = — — (Hidi 3 + fi2<T23 + H.c.), (25) 

and the atom is subject to two quantum jump operators C\ = 1 / 71(731 and C 2 = 

As shown in figure [5] the amplitude-intensity correlation function on the weak 
transition |1) f-)- |3) in a U-system is asymmetric. For r > 0, we observe damped 
dipole moment oscillations, with an additional modulation that we ascribe to the rapid 
dynamics on the stronger atomic transition. For r < 0, the oscillations are more regular, 
however they are damped towards a large almost constant value, and only on a much 
longer time scale (see insert), the correlation function converges to unity. We can account 
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Figure 5. Amplitude-intensity correlation on the weak probe transition of a single 
three-level V atom by invoking the quantum regression theorem (blue solid line) and 
past quantum state (red dashed line). The inset indicates a long (negative) time 
trace approaching unity on the “Zeno time scale”. The parameters are: fR = O.I72, 
fR = 572, 71 = O.OI72, and <p = 7t/2. 


for both of these observations by the backward evolution of the matrix E(t): First, we 
note that E(t) has the excited state 11) (11 as its final state. This state is resonantly 
driven towards the state 1 3), which is, however, strongly perturbed by the laser excitation 
of the |2) |3) transition with Rabi frequency fR- This driving causes an AC Stark 

splitting of the state |3) by fR and detunes the weak transition by 0 2 /2, which emerges 
as a generalized Rabi frequency of the weak and oscillatory amplitude transfer. The 
incoherent rate processes between states |2) and |3) damp these off-resonant oscillations 
on the time scale of a few (72) -1 . The matrix E(t) has not yet reached its steady 
state value, and despite the coherent coupling strength fR, it is only for times longer 
than the time scale ( 71 ) _1 , that the correlation function regresses to the uncorrelated 
product. Within the picture of an evolving matrix E(t), we find this suppression of 
coherent dynamics similar to the quantum Zeno dynamics [23] . observed for evolving 
density matrices m ■ 

6. Conclusion 

Temporal correlation functions witness the non-classical character of light emitted from 
quantum sources, and they have played an immense role in the analysis and in the 
application of quantum states of light. In this article we have revisited the calculation of 
such correlation functions in terms of generalized measurements, and we have discussed 
how the correlations can be interpreted in terms of measurement back action and 
transient evolution of the emitter system quantum state during continuous probing. 

Transients after discrete photon counting events have been well understood as the 
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conventional state dynamics from the non-steady state prepared by the measurement 
back action of the counting event. To account for transients before such discrete 
measurement events, we employed the theory of past quantum states which allows 
calculation of measurement outcome probabilities, conditioned on the knowledge of both 
earlier and later measurements on the system. The practical calculations associated 
with this method are equivalent to the ones of the conventional master equation and 
the quantum regression theorem, but they deal explicitly with the transient evolution 
from a certain initial (final) state, according to an equation with enough similarities 
with the conventional master equation to offer insight into the dynamics. 

We applied the method to amplitude-intensity correlations in the fluorescence from 
laser driven two- and three-level atoms. We recall, that the pertaining correlation 
functions have already been calculated in previous works, but without the interpretations 
offered by the present analysis. 
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